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jy-^ ■ We develop a new formalism for the treatment of gravitational backreaction in the cosmo- 

logical setting. The approach is inspired by projective techniques in non-equilibrium statistical 
(j ' mechanics. We employ group-averaging with respect to the action of the isotropy group of homo- 

geneous and isotropic spacetimes (rather than spatial averaging), in order to define effective FRW 
variables for a generic spacetime. Using the Hamiltonian formalism for gravitating perfect fluids, 
we obtain a set of equations for the evolution of the effective variables; these equations incorpo- 
rate the effects of backreaction by the inhomogeneities. Specializing to dust-filled spacetimes, 
we find regimes that lead to a closed set of backreaction equations, which we solve for small 
inhomogeneities. We then study the case of large inhomogeneities in relation to the proposal that 
backreaction can lead to accelerated expansion. In particular, we identify regions of the grav- 
itational state space that correspond to effective cosmic acceleration. Necessary conditions are 
(i) a strong expansion of the congruences corresponding to comoving observers, and (ii) a large 



negative value of a dissipation variable that appears in the effective equations (i.e, an effective 
"anti-dissipation"). 



1 Introduction 
1.1 Preamble 

The fundamental postulate of modern cosmology is the assumption of a homogeneous and isotropic 
universe. The spacetime must then possess a six-dimensional group of spacelike isometries, i.e., it 
must be of the Friedmann-Robertsn-Walker (FRW) type. However, isotropy and homogeneity refer 
to a coarse-grained level of description: there is significant inhomogeneity at short length- scales. 

Since homogeneity is approximate, one may inquire how inhomogeneities affect the evolution 
of the FRW variables. This question is of foundational interest, as it touches upon the domain of 
validity of the fundamental assumption of modern cosmology HI 13- Moreover, both in early universe 
cosmology (inflation in particular) and in present-epoch cosmology, gravitational backreaction effects 
may play a significant role in the evolution of the universe. In particular, it has been proposed that 
the backreaction of spatial inhomogeneities may be responsible for the apparent cosmic acceleration, 
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so that one would not have to invoke the existence of dark energy OlSl EJ— see also the reviews [J6l [Vj 
and critique j8]|. 

In this paper we develop a systematic approach for the study of gravitational backreaction with 
an emphasis on the cosmological context. Our starting point is the observation that the issue of grav- 
itational backreaction has many analogues in non-equilibrium statistical mechanics. The description 
of a spacetime by a single quantity (the scale factor), obtained from the "averaging" of a generic 
metric, resembles the description of many-body systems through coarse-grained bulk variables (the 
mean field approximation in particular) flUHOldlKEl. A key problem in the study of gravitational 
backreaction is the consistency of the approximation scheme. Problems of this type are specifically 
addressed by the techniques of non-equilibrium statistical mechanics. A transfer of ideas from this 
field of research for the study of gravitational perturbations could prove highly fruitful. 

A second problem in treatments of gravitational backreaction arises from the issue of gauge in- 
variance. A definition of effective FRW variables involves spatial averaging of the inhomogeneities, 
and such averaging can be implemented covariantly only for spatial scalars. For this reason, studies of 
backreaction often restrict to averaging the Hamiltonian constraint, which is a spatial scalar; the av- 
eraging cannot be implemented at the level of the tensor-field- valued equations of motion. Moreover, 
even if one restricts to scalars, the results depend on the choice of foliation, i.e., on the family of spa- 
tial surfaces upon which one integrates. Lack of gauge invariance may also pose a problem in another 
stage, namely, in the implementation of the dynamics. If the dynamics of gravitational perturbations 
are described through gauge-fixing, then there is the danger that that the study of back-reaction will 
lead to results that depend on the chosen gauge [H. 

The formalism we develop in this paper is gauge-invariant. The main difference from previous 
approaches is that we employ group averaging over the isometry group of the FRW spacetime, rather 
than spatial averaging. Group averaging is defined covariantly for any tensor field, and it reduces 
to spatial averaging for scalar quantities. As a matter of fact, the properties of the group-averaging 
calculus allow all group averages that appear in this paper to be reduced to spatial averages of scalar 
quantities. 

In order to avoid the problems related to gauge-fixing, we work within the Hamiltonian formalism 
for general relativity. In this paper, we also restrict the matter content to perfect fluids. In this case, the 
solution of the constraints can be implemented in a fully geometrical way (i.e., without gauge fixing) 
at the level of the Lagrangian. The reduction entails a solution of the diffeomorphism constraints: this 
corresponds in an implicit "selection" of a class of foliations tied to the perfect fluid. The only "gauge 
choice" that remains is that of a time variable. The formalism allows the derivation of backreaction 
equations for any such choice, as long as it can be made consistently over the system's state space; 
we choose time as measured by observers comoving with the fluid. 

1.2 Backreaction and non-equilibrium statistical mechanics 

A common approach for the study of backreaction effects in a cosmological spacetime employs a 
perturbation expansion around the FRW solution. This involves solving the linearized Einstein equa- 
tions around the classical FRW solution; the perturbations are then employed for the construction 
of an effective stress-energy tensor, which, when inserted into the Einstein equations, provides the 
corrections of the FRW evolution. In general, this procedure suffers from consistency problems. In 
particular, the construction of the effective stress-energy tensor is gauge-dependent-hence the result- 
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ing backreaction equations are also gauge-dependent jH. 

Moreover, when large perturbations are taken into account, the accuracy of such an approximation 
scheme degenerates rapidly with time: the solution of the equations of motion with backreaction 
diverges cumulatively from the FRW solution. This means that the perturbations around the FRW 
equations capture less and less of the physics of the system as time increases. Outside the gravitational 
context, such treatments are known to misrepresent backreaction-induced effects such as dissipation 
and diffusion. 

The consistent treatment of backreaction is a major ingredient in most techniques developed in the 
field of non-equilibrium statistical mechanics. The methodology of such treatments varies according 
to the system under consideration. However, all treatments follow a common pattern, which is ab- 
stractly and compactly described in the language of the so-called projection formalism [fT0l[T3l . This 
pattern can be described as a sequence of three steps. 

The first step is the specification of the level of description, namely, of a set of variables that 
provide a coarse-grained description of the system under consideration. For example, in quantum 
Brownian motion, one studies a selected particle interacting with a heat bath of harmonic oscillators 
(environment). The level of description corresponds to the degrees of freedom of the selected particle. 
In Boltzmann's treatment of the rare gas, the system is a collection of weakly interacting particles, 
and the level of description is defined by a probability density on the phase space of a single particle. 
In general, the level of description corresponds to a subspace of the space of functions F(F) on the 
system's state space T. It is represented by a projective map P on F(T). We shall call the variables 
that lie within the range of P relevant variables and ones that lie on the range of 1 — P non-relevant 
variables. For the cosmological perturbations considered here, the relevant variables correspond to 
homogeneous and isotropic field configurations. 

The second step involves a splitting of the dynamical evolution into components in accordance 
with the chosen level of description. Let L t be the evolution operator on the space of states (the 
propagator of the Liouville equation in a Hamiltonian system). 

• PL t P describes the self-evolution of the relevant variables. 

• (1 — P)L t (l — P) describes the self-evolution of the non-relevant variables. 

• (1 — P)L t P describes the coupling between relevant and non-relevant variables. 

The splitting above allows for the derivation of a set of evolution equation for the relevant vari- 
ables: this contains the evolution terms PL t P for the relevant variables and backreaction terms that 
arise from the coupling (1 — P)L t P of relevant to non-relevant variables. The backreaction terms de- 
pend on the state of the non-relevant variables and on their self-evolution in terms of ( 1 — P) L t ( 1 — P) . 
In general, the set of evolution equations for the relevant variables is not autonomous, or not-closed, 
and bears an explicit dependence on the initial state of the non-relevant variables. 

The third step is the derivation of a closed set of equations (e.g., a Fokker- Planck-type equation 
for quantum Brownian motion or Boltzmann's equation for the rare gas). To this end, one introduces 
additional assumptions about the state and evolution of the non-relevant variables. For example, one 
may assume that the state of the non-relevant variables is not significantly affected by the evolution of 
the relevant variables (Born approximation), or that the memory effects in the evolution are negligible 
(Markov approximation), or that the non-relevant variables are "fast" in relation to the relevant ones, 
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or that the correlations between relevant variables are insignificant (e.g., in the truncation of the 
Bogolubov-Born-Kirkwood-Green-Yvon hierarchy in kinetic theory), and so on. Such assumptions 
are necessary for the closure of the set of evolution equations and they often involve the introduction 
of semi-phenomenological parameters that describe the properties of the irrelevant variables. 

To summarize, the general procedure for the consistent construction of backreaction equations 
involves: 

1 . a specification of the level of description. 

2. a splitting of the dynamics and the construction of equations for the relevant variables 

3. additional assumptions about the state of the irrelevant variables that allow for the closure of 
the system of effective equations. 

1.3 Our approach 

In this paper, we apply the reasoning above to the treatment of the backreaction from cosmological 
inhomogeneities. The rationale is that the FRW variables provide a coarse-grained level of description 
for a cosmological spacetime, which is, in a sense, analogous to the description of statistical systems 
in terms of mean-field theory. The backreaction of inhomogeneities is then conceptually similar to the 
incorporation of the effects of the second- and higher-order correlation functions into the mean-field 
evolution Q 

Our primary aim in this paper is to set the basis of a general procedure for the treatment of in- 
homogeneities. The projective formalism can, in principle, be applied to spacetimes with arbitrary 
matter content. Specializing to the case of a gravitating perfect fluid allows us to solve the diffeomor- 
phism constraints in a gauge-invariant way. We then analyze in detail the special case of backreaction 
in a dust-filled spacetime. 

We do not assume a specific form for the "true" spacetime metric, which is to be approximated 
by an FRW spacetime. Rather, we construct backreaction equations in terms of a small number 
of variables that are defined for any spatial geometry, i.e., the backreaction variables are functions 
defined over the full gravitational state space. Hence, the formalism can, in principle, accommodate 
any choice for the "true" spacetime metric. Each choice corresponds to different a region of the 
gravitational phase space, in which the backreaction variables take different values and may generate 
qualitatively different evolutions. 

The method developed here implements steps 1 and 2 described in Sec. 1.2, in a gauge-invariant 
way. The level of description is the subspace of the gravitational phase space with initial data invariant 
under the action of the six-dimensional isometry group of the FRW spacetimes, and the projector P 
is constructed by the group averaging of observables. Evolution equations for the relevant variables 
are obtained through the Hamiltonian formalism. 

The third step involves assumptions about the nature of cosmological perturbations, i.e., it requires 
the specification of a region on the gravitational state space. This ought to be an observational, rather 
than a theoretical input, because the intuitions from non-equilibrium statistical mechanics are not 

'We also note that other techniques from non-equilibrium statistical mechanics (mainly functional methods for quan- 
tum fields) have been employed in the context of early Universe cosmology-see ifTSl and references therein. 
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directly relevant to the cosmological context. For small perturbations, we find generic regimes that 
lead to a closed set of backreaction equations, which can solve explicitly. 

For large perturbations, additional variables appear in the effective evolution equations. It is, 
therefore, more difficult to obtain a closed set of backreaction equations. For this reason, in this 
paper, we only study the kinematics of backreaction in the large perturbations regime. In particular, 
we examine whether large backreaction effects can lead to an effective accelerated expansion. We 
find that there exists a plausible regime in the gravitational state space that manifests effective cosmic 
acceleration. This regime corresponds to initial data (present era) in the gravitational state space that 
satisfy a small number of conditions: these conditions are rather restrictive but they are generic, in the 
sense that they do not require a "fine-tuning" of parameters. We find that acceleration necessitates a 
strong expansion of the congruences corresponding to comoving observers (i.e., an "intrinsic" expan- 
sion of the inhomogeneous regions, in addition to the Hubble expansion), and a large negative value 
of a "dissipation" variable that appears in the effective equations. 

However, the physical relevance of this regime remains an open issue: it is necessary to demon- 
strate that the evolution of initial data in this regime correspond to a cosmological histories compatible 
with observations. Such a demonstration requires a full dynamical treatment for large perturbations, 
and it will be taken up in another work. 

It is important to emphasize that in our approach the relevant variables are determined through 
integration over the spacetime group of isometries. In group-averaging the whole set of points of 
a Cauchy surface of a cosmological spacetime is involved; specific subsets cannot be isolated in 
a coordinate-independent way. For this reason, the method, in its present form, cannot provide a 
definition of an "average geometry" for a generic spatial region. The scope is, therefore, different 
from that of the approaches reviewed in 0. In particular, we do not aim to provide an answer to 
questions such as "how does the Universe look at different scales?" [0Q. 

1.4 Structure of the paper 

In Sec. 2, we introduce the notion of group averaging in the context of general relativity, we study 
its properties and develop the calculational tools needed in the remaining of the paper. In Sec. 3 
we present the Hamiltonian formalism of perfect fluids (in the Lagrangian rather than the Eulerian 
picture), following the treatment of Ref. [fT5l : we show that the diffeomorphism constraints can be 
implemented without gauge-fixing. 

In Sec. 4, we define the relevant variables for the treatment of backreaction, we elaborate on the 
Hamiltonian equations of motion, and we construct the backreaction equations. In Sec. 5, we study 
the evolution of the backreaction parameters in the regime of small perturbations. We specialize 
to dust-filled spacetimes and we identify regimes that lead to closed sets of backreaction equations. 
Large perturbations are taken up in Sec. 6: we find the corresponding backreaction equations, when 
spatial curvature effects are negligible, and we show that there are regimes that correspond to acceler- 
ated expansion. In the final section, we summarize our results, and briefly discuss possible extensions. 
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2 Group averaging 



As we explained in the introduction, our aim is to incorporate the backreaction of the inhomogeneities 
into the evolution of homogeneous and isotropic FRW metrics. To this end, we must construct a 
map that takes generic inhomogeneous tensorial variables (for example, a Riemannian metric) into 
variables of the same type compatible with homogeneity and isotropy. A naive way to proceed would 
be to integrate the tensor field over a Cauchy surface and divide by the volume. However, coordinate- 
invariant integration can only be defined for scalar fields: spatial integration of a generic tensor field 
is not an invariant procedure. In particular, it is not possible to define an average metric through this 
method. 

The method we develop here is based on group-averaging rather than spatial averaging. The 
method is suggested by the fact that the defining feature of the FRW metric is the existence of a 
six-dimensional group of isometries. 

2.1 The basic construction 

Let G be a compact Lie group acting on a compact three-manifold E. This means that there exists a 
smooth map / : G -> Diff(E), such that f 9l o f 92 = f 9l92 , for all g 1: g 2 e G. 

For any tensor field ^4 n "' ln j 1 ...j m on E, we define the group-averaged tensor field (^4 n "' ln j 1 ...j m ) as 

<A''---%,.,J(.r) = J ,1,, if ; r \y - 1 ,....„ {.>■). (i) 

where dfj,[g] is the left-invariant Haar measure on G (normalized to unity). 
The invariance of d/j,[g] implies that 

[/;<^ 1 -"S 1 ...;J](*) = -",,.,,„)(•'•)• (2) 

It follows that there exist vector fields Xa on E that correspond to elements A of the Lie algebra of 
G, such that 

£x A (A h - i " jl ... j J = 0. (3) 

Let T be the state space of Hamiltonian initial data (/i^ , Tc' lj ) for the gravitational field on a space- 
time with Cauchy surfaces of topology E, where hij is a Riemannian metric on E and 7r y the conjugate 
momentum (a tensor density). T = T*Riem(Y), where i?iem(E) is the space of Riemannian metrics 
on E. 

A tensorial variable A that is a functional of and ir tj corresponds to a family of functions on 
the state space T. The map P : F(T) — »■ F(T), defined as P[A] = (A) is a projector on -F(r) and can 
be used to define the relevant variables with respect to the action of the group GonT. 

We also note that the map IT : T — > T, defined as 

n[(h ij ,^)} = [(h ij ),(ir«)}, (4) 

projects onto the submanifold T of T, which consists of initial data invariant under the action / of 
the group G. If G is the 6-dimensional group characterizing homogeneous and isotropic spacetimes, 
T consists of all constant curvature metrics on E and the corresponding conjugate momenta. 
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2.2 Comments 



1 . The construction above applies for a generic (compact) Lie group G and it does not specifically 
require isotropy and homogeneity. For example, G can be a group corresponding solely to special 
homogeneity (Bianchi models). The only requirement is that the group G has a smooth action on the 
three-manifold S. 

2. The group G was assumed compact. (When G implements the symmetry of homogeneity and 
isotropy, E also must be compact: E = S 3 .) However, the idea can also be applied to non-compact 
groups, provided they possess a left-invariant measure dfj,[g]: Let O n be a sequence of open sub- 
spaces of G with compact support, such that O n _i C O n and U n O n = G. Then, for any tensor field 

A %1 ■" tn j 1 ...j m on E we define 

(^■••S,., J(*) = iim ^ /o/^^K/^' - S,-,.^)- (5) 

If this limit exists, and if it is independent of the choice of the sequence O n , then the group 
averaged tensor field (A n '" l "j 1 ... Jm ) is invariant under the action of the group G. The properties 
presented in Sec. 2.1 then follow. 

In the following, we shall assume that the group G is compact, keeping in mind that with suitable 
conditions the results can also be applied to the non-compact case. 

3. From the physical point of view there is an ambiguity in the construction of the projector above. 
The action of a group G on E is unique at most up to diffeomorphisms. If f g is an action of G on E, 
and F a diffeomorphism (that does not coincide with any of the diffeomorphisms f g ), then FofgoF^ 1 
is a different group action on E. 

Group-averaging is equivariant with respect to the action of the diffeomorphism group DiffiT), 

i.e. 

'■"(■■V- J \n..,,„) = ((/•".4)"-'" j! ..,„.)'. (6) 

where (•)' denotes the average with respect to the group action F o f g o F^ 1 . 

Hence, if T is the submanifold of initial data invariant under the /-action of G on E, F*T 
is the submanifold invariant under the F o f o F" 1 action of G. The effective description of the 
group-averaged quantities admits the group DiffiT) as a gauge symmetry. It is therefore necessary 
to reduce the system by removing the gauge degrees of freedom corresponding to the DiffiT)- 
symmetry. This is equivalent to the selection of a specific group action among the class of diffeomor- 
phic equivalent actions. We shall see in the next section, that the presence of a perfect fluid allows for 
a gauge -invariant reduction of the Diff(M) symmetry. 

2.3 The group averaging calculus 

Let A be a tensor field invariant under the action of the group of isometries, i.e., f*A = A for all 
g E G. Then for any tensor field B 

(A®B) = J dpL[g]r g (A ®B)= J d^[g\r g A ® f* g B = J d^[g]A ® f* g B 

= A® j dn\g\r g B = A®(B). (7) 
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This means that an invariant tensor can be taken out of the group averaging operation. The same holds 
if some indices of A and B are contracted. 

Another important property of group averaging is that in a spacetime characterized by homogene- 
ity and isotropy, the group- averaging of a scalar field equals the spatial average of the field on the 
three-surface S. The proof is the following. 

We consider the case of compact three-surface E and group G. Let G x be the stability group of 
x E E, i.e. the subset of G of all elements g E G, such that f g (x) = x. The quotient G/G x coincides 
with E, so that G forms a fiber bundle over E with fiber G x . Considering a local trivialization of the 
bundle <j) x : G — > E x G x , such that 4> x (g) = (f g (x),g'), where g' E G x . Then the measure dfx(g) 
splits as dfJ,-z(y)d/j,Q x (g'), y = f g (x) E E, where d/zs is the G- invariant measure on E. Hence 

(<f>)(x) = J dii(g)(p[f g l (x)} = J dii s (y)4>(y) (J d^ G }j = cj d^(y)4>(y), (8) 

where c — (J d/icj is a constant (it does not depend on </>). Since for the constant function (j)(x) = 1, 
(4>)(x) = 1, c = l/V, where V is the volume of E with respect to the invariant metric. We have 
therefore shown that the group average of a scalar function equals its spatial average over E with 
respect to the group-invariant measure. 

The only properties of group averaging we will use in this paper are the identity © and group- 
averaging of scalar fields. The reason is that we will only encounter group averages of tensors with 
two indices, which, when contracted with the group-invariant FRW metric hij, lead to group averages 
of scalars. It turns out that these are the only averages that appear explicitly in the backreaction 
equations. 

3 The Hamiltonian description of fluids 

The Hamiltonian formalism is well suited for dealing with the problem of backreaction, because the 
level of description associated with homogeneity and isotropy corresponds to a submanifold of the 
canonical state space of a gravity theory. Moreover, in spacetimes with a perfect fluid, the gauge 
symmetry of spatial diffeomorphisms can be factored out completely, due to the special properties of 
the perfect fluid's Lagrangian. Perfect fluid spacetimes suffice for many cosmological applications. 
Dust-filled spacetimes, in particular, are relevant for examining the issue of backreaction-induced 
cosmic acceleration. 

In this section, we present the Lagrangian and Hamiltonian formalism for relativistic gravitating 
perfect fluids. There exist several different approaches [fT4ll : here, we follow an adaptation of the 
formalism presented in Ref. [fT51 . 

3.1 Perfect fluid Lagrangian 

Thermodynamics. The thermodynamic properties of a fluid are encoded in the internal energy 
(Gibbs) functional e(V, S), which expresses the internal energy per particle e as a function of the 
specific volume V and the specific entropy S. The first law of thermodynamics takes the form 

de = -PdV + TdS, (9) 
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where P is the pressure and T the temperature. 

We consider the special case that the internal energy is a function of the specific volume V only- 
i.e., we ignore thermal effects. The internal energy can be written as e(l/n), where n — l/V is the 
number density. The energy density is p = e/V = ne. We note that for an equation of state P = wp, 
the Gibbs functional is 

e = cn w , (10) 

where c a constant. In particular, for w = (dust), we obtain e = c, where c has dimension of mass 
and hence p = cn. 

The matter space. Let Z be the matter space, i.e., a three-dimensional manifold, whose points 
correspond to material particles. These particles are distinguishable, in the sense that each point of 
Z corresponds to a particle of definite identity. The configuration of the fluid is fully determined, if, 
for every point X of the spacetime M, one specifies a particle z G Z(X), whose worldline passes 
through X. Hence, a configuration of the fluid is represented by an on-to mapping ( : M — > Z. Given 
a coordinate system in Z, the mapping is described by three functions 0(^0- 

Globally, if M = £ x R, where £ is a three-manifold, Z must be diffeomorphic to E. Z should be 
a homogeneous space, namely it should carry the transitive action of a Lie group. Z is also equipped 
with a volume three-form v 

v = u(z) dx l A dx 2 A dx 3 , (11) 

which measures the number of particles n(D) within any given region D in Z: n(D) = J D v. 

The pullback (*u is a three-form on the spacetime M. It is closed (d((*u) = 0), since dv = on 
the three-dimensional manifold Z. This implies that the corresponding vector density 

f = - v e^d^C^C 3 , (12) 

satisfies the conservation equation 

d„f = 0. (13) 

The vector density j M is the conserved particle-number current. It can be split into a normalized 
velocity vector field and a particle number density n as 

f = V=gnu», (14) 

whence 



n = -^J-j^fg^. (15) 
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First-order Lagrangian. Having expressed the particle density as a function of the first derivatives 
of the fields z l (X), we write a first-order Lagrangian that reproduces the fluid's equations of motion 



-gp[n\ 



-gne[l/n}. 



(16) 



The three variational equations of L with respect to ( l together with the equation (|I31) lead to the 
conservation equation for the stress-energy tensor |fT5l 



0. 



where 

= -^-g\ pu ^u v + P{g^ + uV)]. 
When the gravitational field is included, the total action 

%,C] = Jd 4 X^(±R-p[nf) 
leads to the Einstein equations G^ v = -T^, where k = 167rG. 



(17) 



(18) 



(19) 



3.2 3+1 splitting 

We next perform a change of variables in the action (fT9l ). This change of variables corresponds to 
a 3+1 splitting of spacetime with respect to the matter space Z. Let t(X) be a time function on the 
spacetime (M,g), i.e., a scalar function on M, such that the condition t(X) = t defines Cauchy 
surfaces S t , and U t E t = M. 

We define a diffeomorphism from M to Z x R as X — > [((X),t(X)]. The inverse map £ : Z x 
R — ► M, (af , t) — > X M = £ M (X) is a functional of C and allows for the standard 3+1 decomposition 
of the spacetime metric g^ u . 

Using the vector fields 

nx) = ^[*po,cpo] (20) 

= ^-[W.CWI (2D 

we express the Riemannian three-metric hij on S t (obtained from the pull-back of g^ u to S t under 

£) as 

hij(t,x) = £f (t,x)£j(t,x)g flu [£ (t,x)]. (22) 

Let n M be the unit normal on the surfaces E t ; then = 0. The lapse function iV and the shift 
vector N l are defined as 

t» = Nn» + £?N\ (23) 
The extrinsic curvature tensor on S t equals 

Ktj := £?£ v 3 V^n u = ^(hij - 3 ViNj - 3 VjNi), (24) 
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with 3 Vj the covariant derivative on Z compatible with the three-metric /^— the dot denotes derivative 
with respect to t. We also note that y/^g = N\fh, where h = dethij. 
The current density of Eq. (fT2l) takes the form 

f = up. (25) 

Substituting Eq. (T25T) into Eq. (fl4]>. we obtain 



'1 - haN'Ni 

" = — 7f, — ' (26) 



where iV* = JSP/N. 

The action (fT9l then equals (up to boundary terms) 

Vh 



S = J dtj d 3 zN^— (KijIC* -K 2 + 3 R- Kne[n}) , (27) 
where K = K^h 1 ^ and 3 R is the Ricci scalar on Z associated to the three-metric h^. 



The action (119]) is a functional on a space Q spanned by the 13 independent field variables g^ v , (\ 
However, after the redefinition of the variables involved in the 3+1 splitting, the action takes the form 
(|27T) that depends on the 10 variables hij, N\ N. The change of variables allowed for the separation 
of the physical from the gauge degrees of freedom. The action (|27l is defined on the space Q re d 
spanned by the field variables (h^, N\ N), 

This change of variables is physically equivalent to the imposition of the comoving coordinate 
conditions. However, no gauge-fixing was involved. The reduction took place at the Lagrangian 
level: the variables h^, N, N l are junctionals of the original variables g^ v , (\ They are coordinates 
on the space Q, invariant, by construction, under the 'gauge symmetry' of spatial diffeomorphisms. 
It follows that (at least locally) Q re d can be identified with the quotient space Q/Dif f(E). 

The construction above does not apply to generic matter Lagrangians: the perfect fluid Lagrangian 
is special in that it allows a gauge-invariant way to select a spatial coordinate system tied on the matter 
degrees of freedom. 



3.3 The Hamiltonian description 

We next perform the Legendre transform of the action (127T ). The momenta conjugate to h^ are 



71 



>J 



—(IC'-Kh* 3 ) 

K 



(28) 



while the variables conjugate to N and iV, vanish identically (they correspond to primary constraints). 
The Legendre transform yields the Hamiltonian 



H= d 3 zN H + fPHi + Vh 



P 



1 - //, ; .Y'.Y< 



x/h 



(29) 
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where 



n { = -2VjTT ij . 

From now on we drop the prefix 3 from the spatial covariant derivative. 
Variation with respect to N l yields 



K 



= up . 



which can be solved for N l : 



N*N 4 



y/ V *F*+WHi 

where F{x) is a function of x = WHij v 2 defined implicitly by the algebraic equation 

Fix) = p [ — = 

\Vh^/l +x/F{x) 2 i 

Substituting Eq. (1331) back to the Hamiltonian we obtain 

U,W -I u/y/h 



H 



d 3 zN 



n + 



P 



Variation with respect to the lapse N yields the constraint equation 



n + 



+ Vh, 



P 



1 + H i H l jv 2 F 2 



0. 



(30) 
(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



The gravitating perfect fluid with the variables above is then a first-class constrained system, 
with the Hamiltonian vanishing on the constraint surface. We see that the use of coordinates tied on 
the matter state have led effectively to an automatic implementation of the spatial diffeomorphism 
constraints. Active diffeomorphisms are not gauge symmetries of the system described by Eqns. (1351 
and®. 

For a dust-filled spacetime (p = cri), we obtain F(x) = c and the Hamiltonian simplifies: 



H 



d 3 zN 



H + Jp 2 + WJU 



(37) 



where p = cv. 
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4 The effective description 



4.1 The symmetry surface 

Having obtained the Hamiltonian for a gravitating perfect fluid, we proceed to the description of the 
effective dynamics in terms of a homogeneous and isotropic metric. 



The Hamiltonian (|35l) is defined on the space T = T* Riem(Z) spanned by the variables h 



Kl 



X 



and n l3 (x). T carries a symplectic form 

Q= I d 3 xd7T ij (x) A6hij(x). (38) 

For concreteness, we set Z = S 3 . Let °hij be the metric of the unit sphere 

°h ij dx i dx j = dx 2 + sin 2 X (d9 2 + sin 2 Odcf). (39) 

This metric is our standard for isotropy and homogeneity in the matter space. It has six linearly 
independent Killing vectors. These vectors generate an action / of the group G on Z. This action 
extends to an action on the space Riem(Z) of three-metrics on Z and on the phase space T. We will 
employ this group action in order to define group averaging for tensors on Z as described in section 
2. 

Let r be the submanifold of T invariant under the action /; we will refer to r as the symmetry 
surface. It consists of pairs (hy,^) of the form h{j = a 2 °hij and 7f y = .p- o °hij. T is 
a symplectic submanifold of T, in which the variables a and p define a symplectic chart, i.e., the 
restriction of on r is 

Qr = dp Ada. (40) 

We also assume that the number density three-form v is invariant under the group action on Z. 
This implies that v = uoV^h, where u is a constant. 

The submanifold T defines the level of description in terms of homogeneous and isotropic vari- 
ables. We shall project the dynamics of the gravitating fluid onto r . The group action / allows us 
to define group averaging for the variables foy, 7T W ' through Eq. (0Q), and thus to construct a projector 
II : T -> T , such that 

n(^,7r«) = «/Ki),<7r«)). (41) 
Hence, the coordinates of any point of T split uniquely as 

(% = (hj + dhij, ^ + 5-K ij ), (42) 
where hj = (hij),^ = and (5^) = 0, (6^) = 0. 



4.2 Choice of time variable 

In order to write the evolution equations corresponding to the Hamiltonian (1351) and the constraint 
(l36l) . one needs to choose a time variable, i.e., to specify the time function t(X) employed in the 
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3+1 decomposition of the Lagrangian — see Sec. 3.2. To do so we must identify a family of timelike 
curves spanning the spacetime M, such that t coincides with their proper time. The equations of 
motion obtained from the Hamiltonian (1351) hold for any choice of time, as long as we take into 
account that the definition of the canonical data depends explicitly in the choice of the time function^. 

The choice of time variable must be made consistently over the whole gravitational phase space. 
This is made possible by the presence of a perfect fluid, because there is a natural choice of time 
associated to comoving observers. Comoving observers correspond to the integral curves of the fluid's 
four-velocity u M — see Eq. (fT4"l) . In matter coordinates, Eq. (fl4)) implies that 

*=_i (»y. (43, 

V-Mm 1 _ n. N . \dtj 

The parameter t corresponds to the proper time along the comoving observers only if 

N = l==. (44) 

1 - N*Ni 

The decomposition ([231) of i M implies that N l = u\ where u l is the spatial component of the three- 
velocity for the comoving observers; u % corresponds to the fluid-particles' deviation velocity from the 
uniform cosmological expansion. N % coincides with Vi = u 1 / \/I + u k Uk, namely, the three-velocity 
of observers along worldlines normal to the hypersurfaces E t . 

Eq. (|33l implies that 

U% = vF(H k H k /v*y (45) 

We note that for a FRW spacetime W = 0; hence, it, = 0: on a homogeneous and isotropic 
spacetime, the worldlines of the comoving observers are normal to the homogeneous and isotropic 
hypersurfaces. Moreover, they are geodesies (since N = 1). 

We next consider a spacetime that is not homogeneous and isotropic. The non-vanishing of the 
vector field it, measures the degree of deviation of the comoving observers from being normal to the 
homogeneous and isotropic hypersurfaces and from being geodesies. 

Test particles move on geodesies of the spacetime metric. This follows under the assumption 
that the particle's backreaction to the geometry is insignificant. However, when we consider a large 
number of particles, we have to take their mutual interaction into account. This results into 'forces': 
the motion of an individual particle is no more a free-fall, and as such it does not correspond to a 
geodesic of the metric. 

When we model an inhomogeneous spacetime by an FRW solution, we essentially assume that 
the world-lines of the particles are geodesies. Hence, we ignore the 'inter-particle interaction' and 
we assume that individual particles evolve under the average fields generated by the whole matter 
content of the universe. In this sense, the FRW description is conceptually similar to the mean field 



2 This means that the time function t(X) must be a functional of the configuration space variables g M „ and £\ Hence, 
the foliation defined from t(X) and ((X) is a functional of the configuration variables. Such "foliation functionals" 
have been introduced in Ref. Ifl6l . The specific functional corresponding to time measured by comoving observers is 
equivariant with respect to spacetime diffeomorhisms. 
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approximation in statistical mechanics. To move beyond the mean field approximation, it is neces- 
sary to consider fluctuations. In many-body systems, the effect of fluctuations is substantial if they 
exhibit strong (long-range) correlation. In analogy, the incorporation of back-reaction in the cosmo- 
logical setting might result in a substantial change over the FRW predictions if there are persistent 
correlations between the properties of the inhomogeneous regions. 

4.3 Equations of motion 

Having chosen the time variable t to be the one measured by comoving observers, we write now the 
equations of motion following by the Hamiltonian (|29l) . We recall that with this choice N l = u l = 

H l /vF(H k H k /is 2 ) and N = ^1 + u* Ui . 

hij = —^(ir ij --nh ij ) + C u h ij (46) 

tJ = _NVh, Rij _ I^i)_2iVK ( « j_ 1 7r7r ii)_Z^^i( 7r ^ + -7r 2 ) 
k 2 ; y/h K 2 ; 2Vh K kl 2 ' 

+ v^W'iV - h ij V k V k N) + C u 7r ij 

+ ^NVh{h ij P(n) + wV[p(n) - P{n))}, (47) 

where P = p — is the pressure, n = a/1 + u k Uk and the Lie derivatives C u read explicitly 

£ u hij = ViUj + VjUi (48) 

C u 7T ij = 7T ik V k U j + 7T kj V k U i -U k V k 7T i3 -7T i] V k U k } (49) 

The initial data satisfy the constraint equation 

$ = £(tt^ - \u 2 ) --R+ [u k u k (p(n) - P(n)) + p(n)\ = 0. (50) 
Note that $ is a scalar on E. 

Symmetric solutions The symmetry surface T is invariant under the Hamiltonian evolution. Set- 
ting hij = a 2 °hij and = 12 ^ a V^h °h^, we obtain u l = 0,N = 1 and reduced Hamilton 
equations, 

" = - K 24^ (51) 

d ,p s , p 2 2ir 2 9 

(« ./„ „ + + 6n 2 Pa) (52) 



dt a 487r 2 a 3 



KOt 



together with the constraint K, g( F^ 4a + — pa 3 = 0. The usual form of the FRW equations then 
follows 



i\\ 2 _ 2 H 



- + a -* = ~p (53) 
at/ o 

a k , 

- = -T^(P + 3P). (54) 
a 12 
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4.4 Effective equations 



We next construct the effective evolution equations for homogeneous and isotropic metrics and 
conjugate momenta tx % K This is equivalent to a projection of the Hamilton equations from the full 
phase space to the symmetry surface r . Let the system be at point (hij, 7T 1 - 7 ) of the constraint surface 
at time t. The corresponding relevant variables are (hij, Tr lj ) = ({hij), (ir tj )); they define a point of 
the symmetry surface. At time t + 5t, the system lies at (h^ + dthij,^ + St-fc^). The effective 
equation of motion is obtained by projecting (h^ + Sthij, ir^ + Stn^) to r . Hence, 

jhij = (hij) (55) 
j t it ij = (>k ij ), (56) 

where hij and tt^ are given by (l46l) and (l47l) respectively. 

Equations (155H56I) are exact: we made no assumptions about the form of the canonical variables 
h^, In order to derive a meaningful effective description, we must specify a region in the system's 
phase space. This is equivalent to a choice of an approximation scheme, according to which the right- 
hand-side of Eqs. (|55H56I) will be evaluated. 

Here, we employ the following approximations. 



1. Perturbation expansion. We first effect of a perturbation expansion around the symmetry surface. 
We point out that we do not perturb around the solutions of the FRW equations, for in this case it 
would be impossible to obtain consistent backreaction equations. 

The first-order terms in the perturbation expansion vanish. This is easy to see for ultra-local 
terms in the perturbations (for example ones that involve (5hij)), because these vanish by definition. 
There are however non-ultralocal first-order terms, such as VkShij, which do not vanish identically. It 
turns out that the contribution of all such terms reduces to integrals over a total divergence and hence 
vanishes for compact spacetimes. 

(This property follows from the high degree of symmetry of the FRW cosmologies. If, for ex- 
ample, we had considered expansion around the symmetry surface for a three-dimensional isometry 
group, the first-order contributions would be non-vanishing.) 

Hence, the dominant terms in the back-reaction are second order to the perturbations. We note 
that to second order iV ~ 1 + ^h^UiUj, i.e., deviation velocities are "non-relativistic". 

2. Ignore ultralocal terms. . The second assumption in our derivation of the backreaction equa- 
tions is that the fluctuations Shij(z) and 5ti^(z) are small with respect to the averages hij(z) and 
7f y (z) in an ultralocal sense. By this we mean that the trace norm of, say, the matrix 8hij(z) must 
be a fraction e << 1 of the trace norm of the matrix hij(z) at almost all points of z. Note that there 
is no ambiguity in the notion of the norm, because hij(x) is defined with reference to the comoving 
co-ordinate system and there is no "gauge" freedom from active diffeomorphisms. However, the per- 
turbations are functions of x: if their characteristic scale is of order / << ljj the action of a derivative 
operator leads to a term of order e/l » e. This implies that the ultralocal terms in the backreaction 

3 The coordinates on Z are dimensionless corresponding to the unit sphere. In terms of distances, the condition I « 1 
translates to fluctuations at a scale much smaller than the scale factor a. 
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equation are much smaller than the contain derivatives of 6hij, 5tx % K We therefore drop all terms that 
involve only ultralocal terms of the perturbations. 

Using the approximations above, we consider the special case of a dust-filled spacetime p = cn, where 
c a constant. We write \i = cv, and /i = cz/ , so that p = p V®h; hence p = c(p/\/h) Vl + UkU k . 
The velocity field u, L satisfies u { = Tti/p. 

We define the tensor field Cj fc through the equation 

W - %V* = Cl k V\ (57) 
where V« is the co variant derivative associated to = (%}. From this definition, we obtain 

C U = \h U (VjShu + V k 5h fl - ViShjj) . (58) 
We shall also use the variables 

A 1 : = h kl C\ l = --^V k (Vhh tk ), (59) 

= (60) 

With the definitions above, W = VjSn^ + ir kl C l kl . Since a term Sn^C^ is of a higher order and 
does not appear in second-order expansion 



W = -2(V^tt ij + Tt ki C l kl ) = -2(V>^ + T77 ^V hX). (61) 



-Mr<i \ _ nix-! £_« i P 

I2n 2 a< 

Expanding Eqs. (|55H56I) to leading order (second) in perturbations around the symmetry surface 
and dropping the sub-dominant ultralocal terms we obtain 

(1 + - ~r« (62) 

1 + ^) + — (1 + ^) 
I not I 





d 

—a = 
dt 


Kp 

24n 2 a 


d 


0- 


np 2 


dt 




48vr 2 a 3 1 




+ 


2n 2 a . 
A ■ 



2/^0— 2 n 
ar 6 a 



u 2 - -T-, (63) 



The quantity A is obtained from the perturbations of the scalar curvature and equals 

A = h ij (Vh(R ij - -h ij R)). (64) 

The quantity u 2 is the mean- square deviation velocity 

u 2 = h lj (uiUj), (65) 

and 

r:=(A*«i), (66) 
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plays the role of a time dependent "dissipation" coefficient. 

We also average over the constraint Eq. (l50l) . Keeping only the second-order non ultralocal terms, 
we obtain 

^_ + ^_ + I«?_^ (1 + Il ? ) = 0, (67) 
9o7r 4 a 4 kot k cr 2 

In Eq. (|67l) 5R is the perturbation of the Ricci scalar 5R = (R) — R, where R is the Ricci scalar 
associated to the metric hij. We find, 



6R=(K i K i -C i ' k C kij ). (68) 

Eqs. (|62l l63l 1671) form a system of equations for the effective evolution on the symmetry surface 
r . These equations depend on four functions of time T(t),u 2 (t), A(t), SR(t), whose form depends 
explicitly on the state of the non-relevant degrees of freedom; they are also functionally dependent 
on the relevant variables a(t) and p(t). Since we have three differential equations determining the 
motion in a two-dimensional surface, one of them can be viewed as a compatibility condition between 
the four functions of time. It is convenient to use Eq. (l63l) for the determination of the independent 
function A(t). Then the effective dynamics on T , can be described by Eqs. (1621 ) and (1671 ). These 
take the form 



where we introduced 



a = f (1 + i« 2 ) - ir a (69) 

<OT = ^ 1 
o ba 2 



247r 2 a; 

In order to study the properties of the backreaction Eqs. (|69l ) and (1701 ), we must find explicit 
forms for the functions T(t), u 2 (t) and 5R(t). We explore different regimes for them in the following 
sections. 



5 The regime of small perturbations 

5.1 Time evolution of the backreaction parameters 

Eqs. (l69l) and (1701) provide the evolution of the relevant (FRW) variables including the backreaction 
terms. They are valid at each moment of time. However, they are not, as yet, dynamical equations. 
The reason is that the coefficients T, v? and 5R that appear in these equations depend on the variables 
a and £. We must find their explicit functional relation, in order to obtain a closed set of dynamical 
equations for the level of description. 

To compute the variables Y = (\ l Ui) and u 2 = (u k u k ), it is necessary to find the evolution 
equations for A* and Uj. In this section, we do so, to leading (second) order in the perturbations, 
ignoring the ultralocal terms. 
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Evolution of u 2 . Using Eqs. (13D) . (1371) and the fact that Ui = Hi/ fM we find that on the constraint 
surface 

«i = -Vi(w*« fc ). (72) 

Eq. (1721) is eracf, i.e., no approximation has been employed in its derivation. We note that the right- 
hand-side in Eq. (|72l) is of higher order to the perturbations than the left-hand- side. Hence, to leading 
order u l is a constant. We do not expect this to be the case in general; however, this result suggests 
that for sufficiently small deviations from the FRW evolution, the deviation velocity vary slowly with 
time. 

To leading-order in the perturbations v? = h l ^(t){uiUj), and since Ui is constant 



Evolution of T. In order to calculate A*, we use Eqs. (l59~l) and (l46l) . We obtain 

-V = -3- + -^Lv k n ki - ^-\7 k (irh kt ) + ^=V k [Vh(V i u k + VV - h ik Viu 1 )}. (74) 
a Vh 2Vh Vh 



Keeping first-order perturbation terms around the symmetry surface and using Eq. (|6T|) . we obtain 

X = -3-y d—j= -^X l - ^(tt/VK) + V k (V l u k + VV - h ki ViU l ). (75) 
a Vh 2Air z a z 2 

Then, to second order in perturbations 

f = - npv 2 + -if) + 392 (76) 

a 2 



where we wrote 



and 



^=(^=V k u k ), (77) 

e 2 ~=lm) 2 -2& lJ ev), (78) 



P = /V« 3 (79) 



is the matter density of the corresponding homogeneous and isotropic solution. 



The tensor i;? - is defined as 



= g + • (80) 

Gjj is the expansion tensor for the (typically non-geodesic) congruence of comoving observers, 
defined with respect to the homogeneous and isotropic metric. 



19 



In linear response theory for non-relativistic fluids, the tensor 0^ is proportional to the fluid's 
pressure tensor Q- This suggests that 6^ incorporates the effects of pressure generated by the "grav- 
itational interaction" of the inhomogeneities. It vanishes for a homogeneous and isotropic solution. 

In the regime that is constant, 6^ is constant, too. It follows that to leading-order in perturba- 
tions 

&{t) = (^^j Wfo). (81) 

Eq. (1761) is a differential equation for T in terms of u 2 and Q 2 , which are known functions of time, 
- and £/a, which is described by the evolution equations and ip, which is, as yet, undetermined. In 
effect, ip does not allow for the existence of a closed system of equations at our level of description. 
The derivation of a closed system of equations requires additional assumptions. 

To this end, we consider different types of perturbations and how the evolution of T differs ac- 
cording to the type. 

a. Dominant conformal perturbations. In this regime, the perturbations of the metric and conju- 
gate momentum are of the form Shij ~ bhij and 5n l: > ~ b'Vhh^, for some scalar functions b and b'. 
Then, the term Vk(nh lk ) in Eq. (1751 ) equals 3Vfc7r fc \ and this leads to the following equation 

f = ^l^r - -Kp^ + 362. (82 ) 
a 4 

Comparison of (l82l) to (1761) shows that for this regime 

i> = ^-f?. (83) 



b. Perturbations not affecting the extrinsic curvature scalar. This regime includes the case of 
transverse-traceless tensor perturbations and transverse vector perturbation for both variables. In 
these cases, 5(ir/\/h) = and hence, ip = 0. Then, 

f = ^l^r - Kptf + 302 (84) 
a 



c. Vector perturbations. In this regime, the terms corresponding to A* and V fc<57r fa dominate in 
the fluctuations of the metric and momentum variables respectively. A two-tensor A %: > corresponds to 
vector perturbations if it can be written in the form 

= E { ^S- + - (/5 • Q)^£) «*(*)> (85) 
q \ y y y / 

where q denotes the wave- vectors, corresponding to eigenvalues of the Laplace equation for the metric 
°hij and (3 l a vector-valued function of q. Assuming that 5tc^ and 5{\fhh % i) are of this form, we find 

^/V^-lff + (86) 
V 2 \ 2 Ka J 
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It then follows that 



V = (u l P tJ ^f), (87) 

where = is the projector to the longitudinal part of the vector fluctuations. 
If the vector perturbations are purely longitudinal, then 

4, = (88) 

At the other end, if the vector perturbations are purely transverse, then we fall back to case b) and 
ip = 0. 

d. Mixed perturbations. In the previous three cases, the perturbations of the metric Shy and 
of the momentum 6n % i were of the same type. In principle, other cases are possible. For example, 
the perturbations Stt^ may be purely conformal and the perturbations Shy longitudinal. This is an 
acceptable region of the gravitational state space, since the only restriction is that the metric hy 
and the momentum 7r lJ satisfy the Hamiltonian constraint. For small perturbations, this regime is 
not dynamically stable (i.e., the system does not spend much time in the corresponding phase space 
region), because the linearized equations of motion couple perturbations of the same type. However, 
this regime might contain be stabilized when large perturbations are considered. 

An interesting case is that of traceless perturbations Sn 1 ^ and longitudinal perturbations Shy, for 
which 

ijj = —(u i P ij \ j ). (89) 



Evolution of SR. The scalar curvature R is a function on the state space, and under Hamilton's 
equations evolves as 



2N 

7T njv " 2 



R = —Ry(n ij - -nh ij ) + tfViR. (90) 



To second-order in the perturbation and ignoring the ultra-local terms, we find 

JR = -^LjR + -^v? + ISRViU 1 ) + —(X% - /cVcf) 

5.2 Autonomous backreaction equations 

In the previous section, we derived the evolution equations for the functions T(t), u 2 (t) and SR(t) 
that enter the backreaction equations (|69l) and (1701 ) for the FRW variables a and £. If these equations 
contained no other variables other than T(t),u 2 (t) and 5R(t), we would obtain a closed dynamical 
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description of backreaction in terms of a small number of variables. We would then have effectively 
substituted the infinite-dimensional dynamical system described by the Hamilton equations (|46|) and 
(1471) with a finite-dimensional one. 

However, the set of equations for the variables T(t), v?(t) and<5-R(t) is not autonomous. It depends 
on histories of additional variables, and these histories can be determined by deriving additional evo- 
lution equations. These in turn will contain additional variables and so on. This is a general feature in 
the derivation of effective equations in systems with a large number of degrees of freedom. A com- 
mon practice in non-equilibrium statistical mechanics is to enforce closure on a system of effective 
equations for relevant variables by considering specific regimes, making simplifying assumptions, or 
introducing phenomenological parameters. 

In ordinary statistical mechanical systems, simplifications occur because we have some knowl- 
edge for the system that allows reasonable estimations about the state of the non-relevant variables, 
like, for example, the assumption of local equilibrium or of molecular chaos in the treatment of gases. 
However, in cosmology our present knowledge of the statistical behavior of the inhomogeneities is 
not detailed enough for this purpose. 

A simplification that allows for the derivation of a set of closed backreaction equations is to 
assume that the term involving the perturbation 5R of the scalar curvature in Eq. £70] is negligible in 
comparison to the other backreaction terms. There is no justification for this assumption, other than 
the simplification it entails to the backreaction equations; in particular, it does not follow from any 
observed properties of cosmological perturbations. However, it suffices to require that the contribution 
of 5R to the effective equations is much smaller than only one of the other backreaction terms. This 
terms dominates over the others as it does over 5R. It is not necessary to specify which one of 
these terms is dominant: we keep them all in the solution of the equations. We must remember this 
restriction when the issue of the validity of the approximations involved is to be considered^ 

Setting 5R ~ 0, the backreaction equations under consideration are 



a = £(1 + l -v?) - l -Ta (92) 

e = ^(i+^)-i (93) 

ba 2 

r = ^r-¥& + ^ + m (94) 

where u 2 and 2 are functions of time given by Eqs. (1731) and (f8TI) respectively. 

This set of equations can be made closed only if if) can be written as a function of the other 
variables. We showed in section 5.1 that this is possible for the specific cases a) and b). We shall 
also show that there are physically reasonable approximations that allow closure also for the case of 
vector fluctuations. 

For vector perturbations it is reasonable to consider an interpolation between the purely longitudi- 
nal and the purely transverse cases. The simplest interpolation involves a simple parameter e e [0, 1], 
which takes the value e = for transverse vector perturbations and the value e = 1 for purely longi- 



4 In fact, it is possible to find regimes leading to a closed set of backreaction equation without assuming that SR 
vanishes, but these regimes involve many additional conditions with no clear physical or geometric meaning. 
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tudinal perturbations. The parameter e is introduced by postulating an approximation 

(u'P^^e^^), (95) 

from which we obtain 

= e -f?. (96) 

In order to understand the meaning of the approximation (|951) and to provide an interpretation for 
the parameter e, we next consider a simple model for the perturbations. 

5.2.1 A model for perturbations. 

Let us assume that perturbations 8hij(x) and 6tt^(x) are concentrated in N non-overlapping regions 
U n of compact support on Z, and that they vanish elsewhere. This implies that u l (x) = Yt,n u \ n )( x )> 
where u) n ^ (x) is a field that has support only in the region U n . We also assume that the regions U n are 
much smaller than the curvature radius of the homogeneous metric hy, so that they are adequately 
described by Cartesian coordinate systems corresponding to the geodesies of fry. We denote by V n 
the volume of each region U n with respect to the metric hij and by v n = V n /V the corresponding 
relative volume (V is the volume of £). 

'Rigid' perturbations. The simplest approximation follows from the assumption that the regions 
U n are 'rigid', in the sense that they can be accurately described by the spatial averages u\ n \ of the 
fields iti n Jx). (The spatial average is taken with respect to the local Cartesian coordinate system.) 

This means that w| n )(x) ^ v\ n )Xu n { x )> where xu„ is a smeared characteristic function of U n . We 
note that this case corresponds to 0^- ~ 0. 
In this approximation 

iV n=l Z 

where 

= V n I d3xpi rtn)( x )xuM = yj (^^l^»(k)i 2 , (98) 
where X(7n(k) is the Fourier transform of xu- 

The value of L u\ n \ depends on the shape of the region U n . For a region with no preferred direction 
(i.e., a ball) Xc/ n (k) is a function only of k 2 and Eq. (|98t yields 

L T, i - la* 

U (n) — g M (n)- 

Then Eq. © leads to Eq. (gU) with e = 1/3. 

A smeared characteristic function corresponds to a boundary that is not sharply defined 
avoid problems arising from differentiation. 



(99) 

. It is convenient in order to 
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A second interesting case is that of fluctuation regions U n with one preferred direction. Let 
be the unit vector corresponding to this direction; xu n is a function of r( n ) ■ k and k 2 . Eq. (1981) yields 

L — i 1 — i 3s n 1 , _ u . nr>\ 

U (n) = ~^— U i(n) + g ( W (n) r M«)K(n)' ( 100 ) 

where s n (x) G [0, 1] is defined as 

In general, s n depends on the relative size of the preferred direction: it takes values close to 1 for 
disk-shaped regions and close to for rod-shaped ones. 

Different assumptions about the distribution of the directions of rL) corresponds to different val- 
ues of e. For example, we may assume that the directions r l r n \ and the shape parameters s n are not 
statistically correlated (i) with each other, (ii) with n y and (iii)along different regions U n . Then, Eq. 
(|97T) yields (|95T ) with e = |. Hence, the shape of the inhomogeneity regions does not affect the value 
of e if its determining parameters can be treated as uncorrelated random variables. 

In presence of correlations, the value of e differs. For example, if the mean velocity of the inho- 
mogeneous regions tends to align along the axis rL^, then L «L) = s n u\ n y Substituting into Eq. (l97l) 
(and assuming no correlations for the value of s n in different regions), we find e = -h Y, n s n > i- e -> e 
can take any value in the interval [0,1], depending on the distribution of shapes for the inhomogeneity 
regions. Similarly, if u 1 ^ tends to lie on the plane normal to r l ,y we find that e ~ 1 — — J2 n 

Expanding (contracting) perturbations. If the mean velocities u % , n \ of the regions U n are small 
compared to the values of ul^ (x), then the regions U n are 'at rest' in the comoving coordinates and the 
field uLJx) primarily contributes to the change of their shape. The longitudinal part corresponds to 
change in volume, while the transverse part corresponds to volume-preserving ('tidal') deformations 
of the regions's shape. The parameter e then corresponds to the fraction of the (non-relativistic) 
kinetic energy of the fluid that corresponds to expansion (or contraction) of the region. In this light, 
the approximation (|95T ) is equivalent to the assumption that the tidal deformations in the regions U n 
are not persistent and they vanish on the mean — hence, that ip receives contributions only from the 
fraction of the velocity field that corresponds to changes of volume in the perturbed regions. 

With similar arguments, we can show that for the case of asymmetric perturbations-Eq. (f89l) . the 
approximation (i^P^A- 7 ) ~ (u l hij\ j ) yields 

ijf = e^r (102) 



The parameter e is therefore introduced as a phenomenological quantity that characterizes the state 
of the perturbations (i.e. the actual three- metric) at a moment of time t. In a sense, it is analogous 
to similar phenomenological quantities that are introduced in non-equilibrium thermodynamics, e.g., 
coefficients of viscosity, heat transfer and so on. The approximation (|95l ) is not expected to be mean- 
ingful over the whole of the gravitational state space, but only in a specific region, corresponding, to 
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the cases considered above, e is expected to change in time, as the perturbations evolve. One then 
needs also assume that such changes occur in a time-scale much larger than the ones relevant to the 
resulting effective equations. 

Either with the introduction of the e parameter, or by restricting to the specific regimes considered 
earlier, Eqs. (I92H94I ) together with Eqs. (1731) and (f8TI) provide a consistent and closed set of backreac- 
tion equations to leading order in perturbations around the symmetry surface. These equations hold if 
the perturbation terms are significantly smaller to the ones of FRW evolution. Applied to cosmology, 
such backreaction terms could perhaps lead to corrections on the order of a few percent from the FRW 
evolution. 



5.3 Solutions 



We next solve the backreaction equations (I92H94I) for the scale factor, assuming that the perturbation 
parameters T, v? and 6 2 are small so that only their first-order contributions to their evolution equa- 
tions are signifcant. We assume that the geometry is approximately flat so that in Eq. (l94l) £ >> 1 
and >> 1; this assumption allows for an analytic solution of the differential equations. 
Using Eqs. (O, (ED and ([94]) we find 



(103) 



a 



a" 



a;' 



= a , u 2 (t ) = Mq and Q 2 (t ) = Q\, for some time t . The parameter / takes the 
| for the cases a, b and c of Sec. 5.2, respectively. 



where we set a(to 
values |, 1 and 1 - 

Using Eq. (1921) and keeping terms of leading order to the perturbations we obtain 



dT 

da 



- 2 -r 

a: 



a 9 / 2 a 7 ' 2 
This is a linear inhomogeneous equation with solution 



a 



a J 



3/2 



1 



Oo 

a 



1/2 



(104) 



(105) 



where h\ 



and bo 



^,and r = r(t 



0) 



3 " ~ ' u v k/j.0 

Substituting into Eq. (1921 and changing into the variable u 



2 



in terms of the constants 



Kfl 



a/a ) 3 / 2 , we find 
1 + Au ~^ + Bu^ + Cu^\ , 



(106) 



.4 



B 



C 



3/2 
«0 



12 



KflQ 



-r + 



2^,3/2 



12- 



6 2 a 



o^o 



(107) 
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Figure 1 : The normalised scale factor a/ao as a function of time t in units of the Hubble time H ; Hq = 
%(t ), for 0q = 0.01, ul = 0.0001, T = -0.1. The dashed line shows the FRW solution. 



We then obtain the solution to the backreaction equations 

k 1 + Au~^ + + C«-V3 " 2V 6 1 0j ' U UJ 

In Fig. 1 a plot is shown of a solution to Eq. (II 101) in comparison to the FRW-solution. The 
two curves diverge at early times. However, in this regime the backreaction parameters are large; for 
examples O 2 grows at small a as ~ a -4 . The approximation employed in the derivation of (II 101) 
is therefore not reliable at very early times and the divergence may be absent in the regime of large 
perturbations. 

The approximation is preserved for later times, and then we see that a(t) essentially coincides with 
the FRW solution. Indeed, let us denote by apRw{t) the FRW solution corresponding to a (to) = «o 
(A = B = C = 0), and by a(t) the solution corresponding to a (to) = olq for non-zero but small 
values of the parameters A, B and C. The difference 5a(t) = a(t) — apRw(t) has the following 
asymptitic behavior 

\8a(t)\ ~ Ca(t). (Ill) 

This implies that the change in the Hubble parameter H (t) = (d/a)(t) is \5H(t)\ ~ C, i.e., the FRW 
solution is stable with respect to small perturbations. 

This result holds for symmetric perturbations (e.g., cases a, b and c of Sec. 5.2). If, however, the 
perturbations are mixed (case d of Sec. 5.2), the conclusions may be different. The solution to the 
backreaction equations for ip given by Eq. (11021) can be obtained through a similar procedure to the 
one leading to Eq. (II 101) : 



L 



(") 3/2 dU - / K/i °(t-t ), (112) 



1 + A'u-^ + B'u- 2 ^ + C'u^r- 2V 6 

in terms of constants A',B' and C that depend on the values of 9o,u and T . In this case, the 
deviation 5a from the FRW evolution behaves asymptotically as 

\8a(t)\ ~ C'[a(t)} 1+e . (113) 
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Hence, for large values of t, the change in the Hubble parameter grows as 8H(t) ~ [«(^)] e - 
Even a small value of e produces an evolution that diverges from the FRW prediction. However, 
an assumption in this derivation is that the value of e at time to remains constant in time, which, 
in the linearized treatment of perturbations, is only plausible for very special conditions at t = t . 
Nonetheless, this result strongly suggests that even small perturbations may make problematic the 
extrapolation of the FRW evolution into the asymptotic future. 

6 Large perturbations 
6.1 The approximation scheme 

Part of the motivation for this work is the idea that the cosmic acceleration, deduced from the super- 
nova data, can be accounted by backreaction effects rather than dark energy. This is only possible if 
the perturbations are large: the corresponding terms must be of the same order of magnitude as the 
ones of FRW evolution. The approximation used in the previous section is therefore not sufficient for 
this case: terms of third and higher order to the perturbations 5hij, 5tt^ from the symmetry surface T 
are expected to become significant. 

The effective equations describing the backreaction of large perturbations depend strongly on the 
region of the gravitational phase space to which they correspond. Inevitably, in order to construct 
such a set of equations for the cosmological setting, it is necessary to have substantial information 
about the nature and behavior of inhomogeneities in the present-day Universe. In absence of such 
information, any model describing backreaction of large perturbations must proceed by more or less 
ad hoc assumptions about the nature of inhomogeneities. Such assumptions are necessary both for the 
implementation of a consistent approximation scheme (i.e., which effective variables can be thought 
of as "large"?) and for the closure of the system of effective equations. 

In this section, we will explore one regime of large perturbations, which corresponds to a rather 
natural generalization of the perturbation scheme developed in the previous section. We do not con- 
struct a closed set of backreaction equations for the general case as in the previous section; the prob- 
lem is substantially more involved. Our study is mainly kinematical: we explore the plausibility of 
the idea that large perturbations around an FRW universe can be held responsible for the observed 
acceleration of the scale factor. 

To simplify our calculations we assume that the deviation velocities are very small, namely that 
v? « 1. Indeed, this condition is born out by observations, if we identify Ui with velocities of 
galaxy clusters — see IITTll and references therein. The term u 2 cannot be, therefore, responsible for the 
substantial divergence from the FRW prediction that would be necessary, in order to account cosmic 
acceleration in terms of backreaction effects. We will drop such terms from our calculations. As in the 
previous section, we shall also ignore the the perturbations of the curvature scalar and consider almost 
flat three-geometries. Consequently, the only large backreaction variable in the effective equations 

is r. 

In order to construct backreaction equations for large perturbations, we employ the following 
approximation scheme. We assume that the perturbations 5hij and Stt 1 ^ are small in the ultra-local 
sense, but that they may be large when acted upon by derivatives. This means that the backreaction 
terms may become large because of the strong spatial variation of the inhomogeneities. In particular, 
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we assume that the perturbations 5hij and 5ir^ are a fraction of order e of the averaged variables hij 
and 7T 1 - 7 (with respect to a matrix norm), and that the variation of the perturbations takes place at a 
scale / on the spatial surface. Expanding to the n — th order of perturbations around the symmetry 
surface, a term that involves m derivatives will be of order e n /l m . The approximation scheme we 
shall use in this section involves keeping all terms such that m — n > 0, which means that e is of the 
same order of magnitude with /. 

6.2 The effective equations 

In the approximation scheme described above, the backreaction equations (1691 - 1701) remain unchanged. 
The higher order backreaction terms are characterized by negative values of m — n. Again, as in the 
previous section, one of these equations can be used to determine the functional dependence of the 
curvature perturbations SR in terms of the other parameters. Dropping the contribution of the v? 
terms, Eqs. (16914701 ) take a simple form 

a = f - ^Ta 

e + 1 = ^o 
ba 

Eq. (1741) for the evolution of X 1 is valid also within this scheme. Together with (1721 ), 
following equation for T 

t = ^^T + 3& 2 ~ + uj (116) 

a 

Compared to Eq. (fToT) . terms u 2 have been dropped (in particular ip ~ 0) and there appears a new 
term. 

u = (V l \ l u k u k ) + 2(V k u i C lik u l ) - (V^ASfc). (1 17) 

In the derivation of Eq. (11 161 ), we dropped ultralocal terms 5hij when they are additive to terms 
hif for example, we approximated (V i\ l UkUih kl ) ~ (V 'iX l UkuJi kl ) . 

To simplify Eq. (11171) . we further assume that the tensor 6^ = |(VjMj + VjWj) is isotropic, i.e., 
that 

(-),, B//„, (118) 

where 6 is a scalar function on Z. This means that expansion dominates over shear in the congruence 
of comoving observers. 

Condition (|1 181) implies that 

uj = -3{eX i Ui). (119) 



(114) 
(115) 

it leads to the 
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The right-hand-side of (11191) corresponds to a spatial integral over Z. Schwarz' inequality then 
applies 

M 2 < 9(e 2 ) (((\ k u k ) 2 )) . (120) 

We define the deviation 5T for X k u k as 

(6T) 2 = ((X k u k ) 2 )-T 2 , (121) 

Eq. (11201) then can be written as 



"I < 3v / e 2 v /r 2 + (ST) 2 ). (122) 



This implies that 



uj = 3ri\/Q 2 TJl + (5T/T) 2 , (123) 



for \r)\ < 1. The parameter r\ is, in general, a function of time. We may view it as a phe- 
nomenological parameter characterizing the present state of fluctuations. It characterizes the strength 
of correlations between the backreaction parameters O 2 and r. If rj = 0, then there is no correlation 
andcj = (Q)(X k u k ) = 0, since (0) = ±(V fc u fc ) = 0. If rj = ±1 then the Schwarz inequality (fT20l) is 
saturated and = aX k u k for some constant a at all points of Z. 

6.3 Cosmological models with backreaction 

With the assumptions state above, the set of equations for the backreaction consists of Eqs. (|1 141 
1 1 1 5 L II 161 11231) and the expressions for ij) provided in the previous section. In general, this does not 
constitute a closed set of equations, because we have not provided an evolution equation for 2 that 
appears in Eq. (11231) . The set of equations only closes for the case of longitudinal 'rigid' perturbations 
considered in section 5.3. However, these expressions suffice for establishing the plausibility that 
cosmic acceleration is possible in the regime of large perturbations we have considered here. 
We assume that the effective FRW spacetime is almost flat, i.e. that Eq. ( II 151 ) implies that 



i = r-^. d24) 

6a 



The following equations then apply 



a 1 \ 2 up 



a 3 / 6 



* = -^-Ir^ + lr 2 -ir, (126) 

a 12 6 a 18 3 



r = -2T- + -r 2 + 3 ( © 2 + r/V0 2 r^i + (st/tY ) , (127) 

a 3 V 

Substituting (11271 ) into (|126| we obtain an equation for the acceleration a: 



a ~~ 12 2a 18 



Q 2 + r]\/Q 2 TJl + (5T/T) 2 



(128) 
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Comparing Eqs. (11251 - 11281) to the standard FRW equations, we note the following. First, the 
constraint equation (|125l) is dissimilar to the FRW equations H 2 = |p irrespective of the theory's 
matter content. Dark energy or a cosmological constant or any other form of matter will be added 
to the energy density term. The backreaction term T is additive to the Hubble parameter H, which 
appears squared in the FRW equations. A negative-valued T implies that the actual value of the 
Hubble parameter is larger than the one predicted from the matter content of the FRW solution. 

Moreover, the presence of the Independent terms does not allow us to bring Eq. (11281) into an 
FRW-like form 

ol k , 

a = ~U if) + 3Peff) (129) 

in terms of some "effective pressure" some P e ff- Again, this indicates how the backreaction terms 
differ in their structure to the ones that can be obtained from any form of matter. The presence of 
6-dependent terms in the right-hand of Eq. (|128l) would suggest that at least this part of the accel- 
eration could be interpreted as some sort of bulk viscosity pressure generated by the "gravitational 
interaction" between the inhomogeneities-see [fT8l for a possible relation of bulk viscocity to cosmic 
acceleration. However, the dependence of the acceleration on pressure is quadratic on 6, instead of 
linear as would be expected from an ordinary viscous fluid. Overall, the structure of the backreaction 
terms are very different from what one would expect from any form of matter. 

We next consider Eq. (11251) at the present moment of time t = t . Defining the Hubble constant 
H = (a/ a) (t ) and Q m = we obtain 

If Vt m < 1, then r(£ ) < 0. We introduce the parameter 7 := —T(t )/(3H ); substituting into 
(11261) we calculate the deceleration parameter q = — 2f (t ) as 

<7o = \ + \lo + To + 3 (^o - ^Wl + ^To/To) 2 ), (13D 
where #0 = V / 6^(t ), and we also defined £70 := —(ST)(t )/(3H ). 

6.4 An exactly solvable regime 

Before proceeding to an examination of Eq. (|131l) . we note that in the regime of "rigid" perturbations 
described in Sec. (5.2), the system of backreaction equations is closed. In this case, 6 2 = 0, and the 
set of backreaction equations (II 141 11151 II 161) reduces to 

a = J- -Ta (132) 

V oat 3 

a 1 o 

r = -2-r + -r 2 . (133) 

a 3 

This system of equations can be integrated to determine the functional dependence between T and 

a. 



30 




Figure 2: Plots of the decelerations parameter qo as a function of Q m , for 77 = 0.95, 0q = 0.4 and different 
values of £70: from top to bottom: £70 = 0.8, 5jo = 1.0, £70 = 1.2. 



r / n x 3/2 

— - = (r/r ) 3 / 2 (a/a ) 3 (r - V^n) (134) 

Eq. (11341) is an algebraic equation that can be implicitly solved to determine T as a function of a. 
Substituting T(a) into Eq. (|132l) then provides a closed differential equation for the scale factor a. 

Since T is negative, Eq. (|133l) implies that it was still negative and with larger negative value in 
the past. Then Eq. (|128l) implies that the acceleration ratio a/a was negative and with larger absolute 
value in the past. It follows that Eqs. (|132l) and (11331) describe a universe that has expanded more 
slowly than the FRW solution. Hence, this regime does not correspond to an evolution compatible 
with the supernova data. 

6.5 Cosmic acceleration 

We return to Eq. (11311) . and examine whether there is a regime in which acceleration is possible 
(g < 0). We note the following 

• For #0 = acceleration is not possible. For acceleration, it is necessary that the ^-dependent 
term is negative-valued, i.e., that 77 > and that < 6 < 7770 v/l + (<57o/7o) 2 - 

• A necessary condition for acceleration (q < 0) is that 77^70 > §; this implies, in particular, that 
if £70 ~ 0, no acceleration is possible. 

Plots of qo as a function of 7 are provided in Fig. 2 for specific values of O , V an ^ ^7o- The 
region of the gravitational state space space region that gives rise to negative q is characterized by 
negative value for the parameter T = (\ k Uk). Moreover, the field 6 = |VjM l must be non-zero and 
positive-valued (by virtue of Eqns. (II 191 ), (11231) and of the fact that 77 > 0). This means that the 
expansion of the congruence of comoving observers must be significant. In other words there must be 
a strong expansion of the inhomogeneous regions, in addition to the Hubble expansion. A substantial 
correlation of the fields and X k Uk is also necessary. 
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The above are strong restrictions, but they do not require any fine-tuning of parameters. They 
correspond to a fairly generic region of the gravitational state space. 

The conclusion above only involved a kinematical study of backreaction. To argue that cosmic 
acceleration is a consequence of backreaction, a dynamical study is necessary. Namely, one must 
demonstrate that the phase space region corresponding to acceleration can be reached from suffi- 
ciently generic (i.e., not fine-tuned) and physically reasonable initial conditions. To do so one needs 
to write a closed set of effective backreaction equations and study the properties of its solutions. This 
will be the topic of future work. 

6.6 Special choices for the background metric 

Many studies of cosmic acceleration as arising from the inhomogeneous nature of the universe have 
considered as a true spacetime various forms of the Lemaitre-Tolman- Bondi (LTB) family of solu- 
tions [3|. These describe a spherically- symmetric cosmological spacetime. In particular, void LTB 
models, namely, models describing a large underdense central region surrounded by a flat-matter 
dominated spacetime provide remarkable agreement with a large fraction of cosmological date, pro- 
vided we are located near the void's center. 

In principle, our method can be straightforwardly applied to LTB-type solutions to Einstein's 
equations. In the canonical formulation, such solutions correspond to three-metrics of the form 

ds 2 = f(x)d X 2 + g(x)(d9 2 + sin 2 6d0 2 ), (135) 

and momenta 

n lJ =diag[p / (x),^ 9 (x),^ 9 (x)/sin 2 ^]sin^, (136) 

where Pf,p g are the conjugate momenta of the variables / and g respectively, and the canonical 
variables satisfy the Hamiltonian contraint in Eq. (|36l) . We assumed that the Cauchy surfaces are 
three-spheres and that the perfect fluid is dust. 

The group-averaging calculus yields allows us to define the variables a and £ of the effective 
description, through equations (hi^h^) = 3a 2 and (n^hij/ V^h) = — 3k£. We find 

a 2 = 7^ r dxifix) sin 2 X + 2g(x)) (137) 
Stt jo 

Z = -T~ r d X [p f (x)+P 9 (x) sin 2 x}. (138) 

OHK Jo 

A study of backreaction for the LTB models would then proceed along the lines described in 
sections 4 and 6. The state space of the LTB model is infinite-dimensional, and a full study of 
backreaction would be beyond the scope of this paper. However, there are some points we can make 
in relation to the approximation schemes we have employed. 

In particular, we have relied on the approximation that the ultra-local terms to the perturbations are 
negligible. This means that the magnitude of the metric perturbations 5 f = f — a 2 , 5g = g — a 2 sin 2 x 
must be much smaller than a 2 and a 2 sin 2 x respectively. The same holds for the perturbations to the 
momenta Spf and 5p g . For LTB solutions that do not satisfy this condition, the approximation scheme 
employed in sections 4 and 6 fails, and a different scheme should be developed. 
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The derivation of Eqs. (11251 - 11281) employed additional assumptions. First, a negligible value of 
the mean scalar curvature (R and of the mean-square deviation velocities (u 2 ). These conditions 
are restrictive, but it is possible to find configurations that satisfy them. However, the large degree 
of symmetry of the LTB solutions does not allow for the isotropy of the tensor 0^. The deviation 
velocity field u % is radial (along the djd\ direction) and depends only on x- As a result the only 
non-zero component of 0y is XX . It follows that Eqs. (1125111281) do not hold in an LTB spacetime, 
because of its high degree of symmetry. 

It is sometimes suggested that the use of a spherically symmetric spacetime as a cosmological 
model can be motivated by the assumption of an implicit averaging of observations over the celes- 
tial sphere. In light of our analysis, it is necessary to point out that such an averaging would most 
probably misrepresent backreaction. An averaging over the celestial sphere could be viewed in our 
formalism as a group averaging over an action of the group G = SO (3). The resulting distribution 
of deviation velocities would correspond to a tensor 0jj, whose only non-zero component would be 
XX as above, which may be substantially different from the "true" tensor 0^, and as such to lead to 
qualitatively different behavior for backreaction. Void models, in particular, need not refer to spheri- 
cal symmetry, and it is to be expected that a non- spherically symmetric void (especially with respect 
to the distribution of deviation velocities) would lead to different predictions from the LTB ones. 

In fact, if we have reasons to believe that a family of solutions to the Einstein equations with a 
high degree of symmetry are good approximations to the "true" spacetime metric, it would be more 
convenient to work directly with these solutions. The method we developed here is intended to be 
used for generic "true" spacetimes: its aim is to identify generic variables that drive backreaction, and, 
subsequently, to relate them to observed quantities. In particular, our method would be suitable for 
dealing with "Swiss-cheese"-like models, which have also been studied in relation to cosmological 
backreaction [5|. In fact, the models for the distribution of the perturbations in Sec. 5.2.1 could 
provide a starting point for such a study, without the assumption of spherical symmetry for the void 
regions. 

Finally, we note that Eq. (11251) leads to a relation between the Hubble factor ~ and the quantity 

z — ao/a — 1, of the form 



where 7(2) = — T/(3H ) is a function of z, such that 7(0) = 1 — \/VL m . In effect 7(2) plays the role 

of a redshift-dependent "anti-dissipation" coefficient. 

In contrast, the relation between the Hubble factor and the redshift for dark energy is of the form 



for some function f(z). It is important to emphasize the structural difference between the two ex- 
pressions for H(z): an "anti-dissipation" term T leads to a luminosity-distance relation that differs 
strongly from the one obtained by most reasonable energy-density terms. For this reason, it is con- 
ceivable that good fits to the supernova data may be provided by simple functional expressions for 
7(2) that do not correspond to accelerated expansion. 

The consideration of perturbations of the scalar curvature would lead to a modified expression of 




(139) 





(140) 



Eq. C[39]) of the form 




(141) 
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where r(z) = —6R/(6Hq) incorporates the effect of curvature perturbations. Unlike the 7(2;) term, 
the curvature perturbations are additive to the energy density for dust; their presence in H(z) is like 
that of a dark energy term. 

However, one should be careful before using the variable z as a measure of the observable red- 
shift. The relation between redshift and the scale factor should also take into account the presence of 
inhomogeneities. Averaging over the inhomogeneities is expected to introduce additional terms in the 
relation between the physical redshift and the scale factor. The derivation of such a relation, through 
group-averaging of a generic spacetime is the necessary next step, before to attempt to relate the dy- 
namical equations derived from the present method to observable quantities. This implies that the 
relation between the Hubble factor and redshift will be of the form (11401) . where z will be a functional 
of the physical redshift. This would seem to imply an even stronger divergence from the predictions 
of dark energy models. 

7 Conclusions 

In this paper, we developed a general procedure for the treatment of backreaction in cosmological 
spacetimes. The key ingredients to the formalism has been the averaging with respect to the isometry 
group of FRW cosmologies. This allowed for the construction of a projective map in the gravitational 
phase space for the gravitating fluid, which was used in order to provide a consistent and gauge 
covariant treatment of gravitational backreaction. For dust-filled spacetimes we identified and solved 
the backreaction equations for small perturbations and we identified regions of phase space, in which 
accelerated expansion is possible if the perturbations are large. A dynamical study of the case of large 
perturbations (aiming to construct explicit solution of backreaction equations) will be undertaken in 
a latter publication. 

Some comments on the potential generalizations of the method are in order. We have exploited 
here the properties of the Lagrangian formalism of perfect fluids, in order to factor out the gauge free- 
dom corresponding to spatial diffeomorphisms. In principle, the method can be applied to spacetimes 
with matter content other than a perfect fluid, by taking into account the group-average of (combina- 
tions of) the diffeomorphism constraint. There are, however, issues related to the gauge-invariance of 
the perturbation expansion around the symmetry surface that need to be explored. 

Here, we described backreaction in terms of an autonomous set of evolution equations, construct- 
ing an effective dynamical system with a small number of variables. An alternative approach is to use 
probabilistic arguments: the effective equations would then involve a degree of stochasticity due to 
our ignorance of the detailed state of the perturbations. Models, such as those of Sec. 5.2.1, may be 
useful in this approach, the behavior of perturbations being encoded in a small number of parameters 
characterizing an effective probability distribution for the properties of localized perturbations. 
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